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Abstract 

After reviewing the general scaling properties of aging systems, we present a numerical study of the slow 
evolution induced in the zeta urn model by a quench from a high temperature to a lower one where a 
condensed equilibrium phase exists. By considering both one-time and two-time quantities we show that 
the features of the model fit into the general framework of aging systems. In particular, its behavior can be 
interpreted in terms of the simultaneous existence of equilibrated and aging degrees with different scaling 
properties. 
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1. Introduction 

Slow relaxation [lj-jlj is the process whereby a system evolves without reaching equilibrium in finite 
times. Physical realizations of this phenomenon are observed in many materials, for instance ferromagnets, 
binary or complex fluids and glasses, when they are subjected to a drastic change of thermodynamic control 
parameters, as in a temperature quench or a pressure crunch. In the large-time regime these systems exhibit 
typical features of equilibrium, such as the absence of macroscopic currents and the constancy of one-time 
observables, together with hallmarks of non-stationarity, like the lack of time-reversal invariance and the 
dependence of two-time quantities on the system’s age t w , i.e., the time elapsed after the initial change of 
parameters. 

These apparently contrasting features can be interpreted in a coherent way if the degrees of freedom 
of the non-equilibrium state can be classified into groups with different properties [7, g}. In the simplest 
cases, including ferromagnetic systems and the model considered in this paper, this amounts to identify 
those degrees whose evolution is sufficiently fast and unrelated to the global non-equilibrium condition of 
the system to rapidly self-equilibrate with respect to the new value of the control parameters. These fast 
degrees, whose collection we will indicate with {ip} in the following, are responsible for the appearance of the 
typical features of equilibrated systems discussed above. The rest of the system - constrained degrees which 
are not able to efficiently equilibrate with respect to the new control parameters - are a slow component, 
denoted here with {a}, which gives rise to the aging properties, namely non-stationarity and lack of time- 
reversal invariance. 

This separation into two classes is fundamental to understand the properties of aging systems: the 
possibility to disentangle global properties into their very different contributions has represented a ground¬ 
breaking tool towards the recognition of the ultimate scaling structure underlying the evolution. However, 
despite the relevance of the subject, a precise operative definition of fast and slow degrees is, in general, 
difficult to adopt and a convincing proof of their existence has been given, as to now, only in few cases [gj]. 
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In this paper we consider a simple paradigmatic example of aging dynamics, the zeta urn model loi- 12|. 
where fast and slow degrees can be easily identified, their different properties can be studied, and a clear 
picture of their mutual interplay can be drawn. The scaling properties of {ip} and {a} are explicitly worked 
out and the relative contribution to the behavior of different observables, either dependent on a single time 
or on two times, is discussed. We also discuss how different and a priori unrelated operative definitions of 
{ip} and {cr} lead to the same overall properties, confirming that the separation of degrees is a robust and 
general property of aging which does not depend on the details of the procedure. 

The paper is organized as follows: In section [2] we review on general grounds the behavior of one-time 
and two-time quantities, showing how they differently encode the non-equilibrium properties. In section [3] 
a general identification of slow and fast components is made, their statistical properties are worked out, 
and their mirroring in physical observables is analyzed. In section [1] we introduce the zeta urn model that 
will be considered in the remaining of the paper. The result of a numerical study of the model is presented 
in the following section [5j where we explicitly show a consistent way to identify the sets {ip} and {a} and 
explore their contribution to one-time and two-time quantities. Section [G] contains our conclusions and the 
discussion of some open issues. 


2. One-time and two-time quantities 

Let us denote with n = rip (t) + n a (t) the number of degrees of freedom which are separated into the 
fast set {ip} and the slow one {cr}. Since the system is approaching equilibrium, where by definition out-of- 
equilibrium modes are absent, the measure of {a} is bound to decrease as the evolution proceeds, 

lim n a (t) = 0. (1) 

t—> OO 

This occurs independently of the capacity of the system to reach equilibrium in a finite time (interrupted 
aging) or not. In the following we will always consider the second case, where equilibration is never achieved, 
which requires the thermodynamic limit n —> oo to be taken before the large-time limit t —> oo. 

The organization of the internal modes into two sets is reflected in the behavior of any physical observable 
O. In the case in which {ip} and {a} can be considered statistically independent, an additive form 

o = Op + O ct (2) 


is found, with obvious notation. 

Let us consider one-time quantities first, as for instance the internal energy, the pressure or the density. 
Due to the equilibrated character of {ip}., Op reaches almost immediately its time-independent equilibrium 
value O eq . For example, taking Oit) as the energy, O eq = E eq represents the internal energy that the system 
would possess in the final equilibrium state. If O is sufficiently regular, in a sense that will be more precise 
at the end of this section, since n a (t) —> 0 for sufficiently long times, one generally has 

lim O a (t ) = 0. (3) 

t—>oo 

For this reason, if t is large enough and O eq is finite, the contribution of 0 0 represents only a small vanishing 
correction and this makes one-time observables basically indistinguishable from their constant equilibrium 
value. This however does not mean that the system is equilibrated, but only that the contribution of the 
<t-s is negligible in the one-time quantities. Clearly, if O eq is known the small aging contribution O a can be 
evidenced by subtracting O eq from O. 

The non-equilibrium aging character is instead more manifest in the behavior of two-time quantities 
(again sufficiently regular, see below) 0(t,t w ), with t>t w , such as correlation functions or responses. In 
the asymptotic time-domain in which the age t w of the sample is large, the typical timescales ip and t a 
over which fast and slow degrees decorrelate are widely separated. Indeed, while t a = t a (t w ) is usually an 
increasing function of t w , fp is generally microscopic and age-independent, again due to the fact that {ip} 
is basically in equilibrium. 
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Figure 1: The behavior of C^,(r), C a (t , t w ) and C(t. t w ), as r = t — t w is varied, is schematically shown on a log-log plot. The 
left panel is for a ferromagnetic system, the right one for the urn model. For each function two waiting times are considered, 
the curves corresponding to the larger t w decay later (the curves for Cb, collapse on a single one because this quantity is 
stationary). Note that, as discussed in the text, for a ferromagnet C G (t,t w ) only depends on t/t w . 


This is made evident by the fact that a two-time quantity 0(t,t w ) can probe either fast or slow degrees 
according to the direction along which the large-time limit is approached in the (t,t w ) plane. Specifically, 
if we let t w —> oo while keeping t = t — t w finite, since r <C t a the slow degrees act adiabatically and 
can be considered as static. In this time sector, denoted as the short time difference (or quasi-equilibrium) 
regime, one sheds light on the dynamical behavior of fast modes. Indeed which depends only on 

t = t — t w due to the stationarity of tp, converges to its asymptotic equilibrium value limi-^oo O^ = O eq 
in a rather short time r ~ in which O a cannot change appreciably and remains fixed to its equal-time 
value O a (t w ,t w ) = q. This is pictorially shown in the left panel of Fig. |TJ where the behavior of a typical 
two-time quantity, the autocorrelation function C(t,t w ) of a spin system (see the definition in Eq. ©), and 
its contributions C^(t, t w ), C a (t, t w ) are sketched. The short time difference sector is on the left part of the 
picture. 

Conversely, letting again t. w —> oo but keeping t/t w (or the ratio £(t)/£(t w ), where £ is a monotonously 
increasing function, in specific systems) finite, one has r ~ t w ~ t a : This means that the analysis is focused 
on the timescale of the slow (or aging) components. In this time sector, usually denoted as aging regime (on 
the right side of the left panel of of Fig. [[]), the fast degrees act as a stationary background and the time 
evolution of the slow modes is probed. Indeed, from r ~ t a onward also O a starts changing and O moves 
from the plateau 0(t,t w ) ~ q to the final value linv^oo 0(t, t w ). 

Generally, a good separation of the degrees can only be operated for very large times. The previous 
discussion shows that, in this case, two-time observables are quantities particularly fitted to highlight the 
role played by if) and a. 

We mentioned before that lim^oo O g = 0 occurs for sufficiently regular observables O. We now precise 
this point. Consider, in order to keep the discussion simple, a quantity that can be expressed as O = )Tb o^ l \ 
where the index i runs over the microscopic degrees of freedom and o 1 ' 1 - 1 is the contribution of one such degree. 
This is usually the case when dealing with extensive quantities. Then we can write O a (t) = Oa\t) = 

n a {t)o a {t), where o a (t) is the average contribution to O g provided by one slow degree. If the time-dependence 
of o a (t) is such that lim^oo (f)o ff (f) = 0, then lim^oc 0(t) = lim t _ > , 00 (t) = O eq and we call O an 

equilibrating quantity. Otherwise, the contribution of O a does not disappear and lim^oo O ^ O eq . 

Depending on the system considered and on the nature of O either situation can be observed. In 
most cases o is weakly dependent on time, and the corresponding observable is an equilibrating one. For 
instance the energy of a non-disordered magnet quenched to a temperature Tf below the critical temperature 
converges to the value E eq (Tf) provided by the standard equilibrium statistical mechanical calculation. 
However there are also some notable exceptions, as for instance the response function %(t) (a quantity 
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deeply related to the magnetic energy) of a disordered one-dimensional magnet quenched to Tf = 0, as 
described by the Ising chain with a random magnetic field of strength h. In the limit h —¥ 0 one finds fl3- 17 
li m t->oo X(t) 7^ Xeq- 


3. Identification of stationary and aging modes 


A possible identification of fast and slow degrees, based on arguments similar to the ones discussed in 
the previous section, can be given along the lines of Ref. [s]. Referring to a system described by a collection 
Si = ±1 of spin variables one defines the slow part of Si in terms of the running time-average 

~ At J t dzS i ( Z ) ( 4 ) 

where 

<C At <C t a . (5) 

The fast part is then obtained by subtraction 

ipi~ Si(t)- (6) 


The idea in Eq. © is that, since At is chosen much larger than the contribution to of the fast modes 
self-averages to zero (in a more general situation the fast degrees may average to a finite constant value) 
and we are left precisely with the contribution of the slow modes. 

The condition At -C t a guarantees that the slow modes are frozen on the timescale At, therefore the 
running average in Eq. 0 does not mix its value at time f with the future evolution. 

Clearly, as it is, Eq. © is tautological since we usually don’t know t a unless the slow degrees are defined. 
However, using the properties of two-time quantities described in Sec. [2] and sketched in the left panel of 
Fig. [ 1 ] the condition 0 can be translated into 


1 

A t 


r t+At 


C(z, t w ) dz ~C(t + At, t w ) 


or equivalently 


for the autocorrelation function 


h 


't w -A t 


C[t , z) dz ~ C(t, t w — At), 


C(t,t w ) = (Si(t)Si(t w )). 


( 7 ) 

( 8 ) 
( 9 ) 


This can be understood looking at the schematic plot on the left panel of Fig. [T] Indeed, let us consider 
the three possible choices of t: i ) t w < t < t w + t^,, ii) t w + t^> < t < t w + t a and iii) t > t w + t a . 
Starting from ii), in this case the correlation function is constant and equal to q (to be precise this occurs if 
t w + < t < t w + t a — At but t w + t a — At ~ t w + t c r because of d5|)) and therefore equations 0 and 0 

are obvious. On the other hand, considering for simplicity Eq. 0 (for Eq. @ one can proceed similarly), 
in case i) C(z, t w ) is not a constant in the interval z £ [t w , t w + 1^\ but this integration interval is negligible 
if the condition At t^ in Eq. 0 is met. Finally, Eq. 0 is true also in the last case iii), because in this 
stage C is decreasing but this occurs on timescales of order t a , while we have chosen At <C t a . 

Letting t w = t in Eq. 0 , with t chosen as in ii) above, one has 


C(t + At, t ) ~ q, (10) 

which very transparently appears as a direct manifestation of the choice 0 . The quantity q in this case 
amounts to the so-called Edwards-Anderson order parameter. In the previous argument and in what follows 
it is understood that the symbol ~ is correct in the large-t^, limit, and the corresponding equalities become 
exact (i.e. with ~ replaced by a strict equality =) for t w —> oo. 
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Let us notice that the modes defined through Eqs. flUt are statistically independent. Indeed any 

cross-correlation between them vanishes 


(c Ti(t)tpi(t w )) = (■)} = 0, 


( 11 ) 


as can be easily proved 


= (ai(t) [Si(t w ) - ai(t w )]} 

i nt-\- At -i pt-\- At pt w 

= mL c ^ d *-w?L m c ^ v)dzdv 

~ C(t + At, t w ) - C(t + At, t w - At) ~ 0, (12) 

where we have used Eqs. (EH) and the last passage can be obtaind by reasoning as below Eq. © enforcing 
the properties of C. The second relation in (fill) , namely ( ipi(t)<Ji(t w )) = 0, can be proved analogously. 

Due to the statistical independence ED, the correlation splits into 

C{t,t w ) = C^(t) + C a (t,t w ), (13) 


where we have assumed that the ipi’s are stationary variables. The schematic behavior of Cy, and C a is 
shown in the left panel of Fig. [L] 

While inherits the properties of the target equilibrium state, the aging part O a displays new features 
characteristic of non-equilibrium, as a remarkable gauge invariance called dynamical scaling. This can be 
expressed as the fact that the flow of time can be fully taken into account by changing the units £ of 
measure of certain dimensional physical quantities. Assuming, as it usually (but not necessarily) is, these 
to be lengths, one has invariant forms by measuring lengths in unit of the time-dependent unit of measure 
£{t). For instance, in ferromagnets the equal-time correlation function G(r,t ) = ( Si(t)Sj(t )) between spins 
i,j at distance r obeys 


Ga{r,f) = g 


ml ’ 


(14) 


where g is a scaling function. The clear meaning of this is that configurations of the same system at different 
times look the same on average if we rescale space by the characteristic unit measure of lengths, which in 
this case has the very transparent meaning of the size of the growing ordered domains of the two symmetry 
related equilibrium phases at T = Tf. More generally, for a one-time quantity which depends on space, 
scaling implies 


O a {r,t) =r\t)d 


mV 


(15) 


Eq. (THU) being a particular case with 5 = 0. 

In a similar way, for the autocorrelation of a ferromagnet one has C a (t,t w ) 
particular case (again with b = 0) of the general scaling 


c 


m 


which is a 


O a (t,t w ) = £ b (t w ) d 


' m - 

_£{tw). 


(16) 


In ferromagnetic systems with a scalar order parameter^ like those described by the Ising model, the two 
classes of degrees {t/’},{CT} have a natural interpretation 18]. After a temperature quench domains of the 
two equilibrium phases with magnetization ±M form, compete and grow. In this case <r,; and ifi can be 
thought as the background magnetization of the domains, which goes from +M to —M passing from zero 
upon moving across a domain’s boundary, and the fast thermal flips of spins on top of this background, 
respectively. 

Interestingly, a separation of degrees of freedom into two classes with well separated time-scales can be 
exhibited in an exactly solvable model describing a ferromagnetic system with a vectorial order-parameter 
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with A f components, in the large -Af limit [9]. The order parameter cf>(x,t) - which can be viewed as a 
coarse-grained spin - can be explicitly decomposed as <j>(x, t ) = ip(x, t ) + <j(x, t) in an exact analytical way, 
and the statistical independence CD is proved. 

The previous considerations show that, in some cases, the splitting of the degrees can be achieved. 
However, in the previous examples the two kind of modes are spatially mixed because both components 
are uniformly spread around the whole system. In the rest of the paper, instead, we will concentrate on a 
simple aging system, the zeta urn model jiol - tlii . which, besides the advantage of being analytically 

tractable, has the remarkable property that {tjj} and {a} are segregated in different parts of the system. 
This allows to concretely visualize them and offers a more intuitive interpretation. 


4. The zeta urn model 


The zeta urn model is defined by a collection r* € {0,1 ,,M} of integer variables defined on a network 
of N nodes i. Usually the nodes are identified as boxes or urns inside which r. t marbles, whose number 
M = r J i is fixed and conserved, can be accommodated. A cost function, or Hamiltonian, is introduced 
as the sum of independent contributions from all the nodes 

N 

H{{ri}) = '52 E i, ( 17 ) 

Z=1 

where 

Ei = ln(r* + 1) (18) 

is the (energy) cost to accommodate marbles inside the i-th box. In a statistical mechanical approach 
one introduces a canonical ensemble at the inverse temperature /3 and configurations {r*} are weighted by 
the Boltzmann measure 

V({r x }) = Z~ 1 (f3, TV, M)e- /3 ' H({r * }) (19) 

where 

Z(p, N,M) = Y j e“ /m({ri} ) (20) 

In} 

is the partition function. An exact computation shows the existence of a density-dependent critical tem¬ 
perature T c = (3~ l . Above such temperature the probability P(r) of any of the identical urns to contain r 
marbles decays exponentially with r, while 

P(r)~(r+1)"*, (21) 


with /3 = /3 C , at the critical point. In the low-temperature phase a macroscopic condensate appears: the 
occupation probability P(r) is still given by Eq. (12T1) . with the current value of /?, for all the boxes except 
one, which contains an extensive number of marbles. 

The equilibrium measure is invariant under a dynamics that respects detailed balance. In the following 
we will consider a kinetic rule where the number of particles is conserved and the transition rate w(r.i,rj) 
to move a particle from a non-empty urn i to another j. chosen randomly in the network and containing 
respectively r, and r 3 particles before the move, is given by the Metropolis expression 

w(ri,Tj) = min [1, exp(— j3AE )], (22) 


where 

A E = ln(r-j) — ln(r 4 + 1) + ln(r^ + 2) — ln(rj + 1) (23) 

is the variation of the energy due to the move. Notice that the present kinetic rule has a mean field nature, 
in the sense that particles can be transferred between any couple of nodes. 

The dynamics of the present model was solved exactly in different conditions in Ref. 22] where many 
quantities were computed. Here, since we are interested in separating degrees of freedom with different 
properties, we study it numerically. 
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5. Numerical study of the model 


We performed a series of simulations of a system of TV = 500 urns, which are initially prepared in an 
infinite temperature equilibrium condition with particles distributed according to 


P(r) 


{P + IY+ 1 ' 


(24) 


where p = M/N is the density which, in the following, we will set to p = 2. From this initial condition we 
quench the system to the final inverse temperature (3f = (/cbT/) - 1 = 6 > /3 0 , meaning that we evolve the 
dynamics with this value of /3 in the transition rates m • Observables are averaged over a non-equilibrium 
ensemble of 3 x 10 5 realizations of the process with different initial conditions and thermal histories. 

The target equilibrium state with a macroscopic number of marbles condensed in a single urn is ap¬ 
proached by a dynamical process that can be roughly divided into three stages. In an early stage, for 
0 < t <C TV (where t is the time elapsed after the quench, measured in units of Monte Carlo sweeps), a 
segregation occurs between urns with a low occupation number, that we will denote as the gas phase in the 
following, and a liquid phase containing boxes occupied by many marbles. For larger times TV <C t <C TV 2 a 
scaling regime sets in during which the gaseous regions remain on average stationary while urns of the liquid 
phase compete and, as an effect of this, their number diminishes. Eventually, for t N 2 a finite-size effect 
sets in, the symmetry between the nodes is spontaneously broken, a finite fraction of marbles is stored in a 
prevailing urn and the rest of them is dispersed among the remaining TV — 1 boxes. 


5.1. Occupation probability 

This whole pattern of behavior is mirrored by the evolution of the occupation probability P(r, t), as it is 
shown in Fig. [2] From the initial distribution, P(r, t) starts to progressively form a hollow at an intermediate 
value r mln of r, signalling that the segregation process is started, and around t ~ 10 2 a secondary peak is 
developed around r = r max = 10. From this time onward the peak moves progressively to larger values of r 
and becomes more pronounced as compared to the hollow minimum. We will say that urns with r < r m i n 
are in the gaseous phase and the remaining ones in the liquid one. 

Around t ~ t sc ~ 1.5 • 10 2 the second regime is entered, where the form of the different curves do not 
seem to change except for a rescaling of both axis. 

Eventually, around t = 10 5 a sharp peak located at the upper right of the spectrum starts growing. This 
is the hallmark of the final stage of the process where the finite size is felt and all the condensed marbles 
fill a single urn. 

In the following we will concentrate on the second regime where the scaling properties are observed. Let 
us stress that, if the thermodynamic limit TV —> oo (with fixed p) is taken from the onset, this stage becomes 
truly asymptotic. As discussed in Secs. [Q [21 03 this is the regime where the two classes of degrees {tjj} and 
{a} coexist. 

We show now that these two sets can be consistently associated to the gaseous and the liquid urns. In 
order to do that, we observe that - as it can be seen in Fig. [2] for r < r m i n the curves are time-independent 
and all collapse on the equilibrium power-law distribution m■ Therefore urns with r <C r m i n are already 
in an equilibrium state and this is the characteristic property of the {ip} component. Conversely - as it 
is shown in the inset of Fig. [2] by plotting tP(r,t ) vs r/t 1 / 2 , for r r min the curves at different times 
collapse on a single master curve. This means that the form 

is obeyed, where [l9| 

£(t) = t 1/2 

is a time-dependent characteristic marble number and p is a scaling function. 

Equation (1^51) is a typical scaling property, formally identical to Eq. (Hsu, with the difference that r is 
not a space coordinate but, say, the height of the stack of marbles in an urn and £(t) the unit to measure it. 


(25) 

(26) 
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Figure 2: The occupation probability P(r,t) is plotted against r + 1 on a log-log plot for (3f = 6 and p = 2. The dashed 
violet line is the law (r + 1) —6 , namely the equilibrium occupation probability J21D of the target equilibrium state. In the inset 
tP(r,t) is plotted against r/t 1 / 2 . 


This suggests the quite natural identification of the highly populated urns with r > r m i n as {cr} (hence the 
subscript in Peril]. Notice that the property (E5ll implies that the number n a {t) ~ °f liquid 

urns decreases as 

n a (t) ~ ({t ) _1 , (27) 

as it is expected for an aging component according to Eq. ©• Since the number of marbles in the condensate 
is roughly constant during the scaling regime, the shrinking of n a (t) entails a correspondent growth of the 
occupation number of a liquid urn proportional to i(t). This justifies the physical interpretation of £(t) as 
the natural unit of measure of r. 

The scaling (E5l) implies also that linp^oo P a (r, t) = 0, that is the property ©, hence P(r, i) is an 
equilibrating quantity. However quantities like the fluctuation ( r 1 2 ) a — (r)„ of the number of liquid marbles 
(where (r) a = P(r,t) and similarly for (r 2 ) CT ) are not equilibrating, as it can easily be checked using 
Eq. (1551 . 

5.2. Energy 

As a different example of a one-time quantity we have considered the energy density e(t) = E(t)/N = 
( H)/N of the system, where the Hamiltonian is given in Eq. (1171) . Following the decomposition fl2j) we write 
e(t) = e^,(t)+e a (t) and, using the non-interacting character (fT71) of the model, we define - g r^i Ei/N 

and e a = J]je{cr} Bi/N (the two sets {?/>} and {a} are identified as discussed above). These quantities are 
shown in Fig. [3] e^(t) converges to the final equilibrium value e eq {Tf) = in a time of order 


1 More precisely, one should say that the highly populated urns contain also, besides the condensate {rx|. a fraction of the 

gas phase. The latter is represented by the few marbles hopping from/to the scarcely populated urns. However, for sufficiently 
low temperatures this gaseous fraction is very small (less then one ball in a condensed urn, in our simulations), so we can 
neglect it and identify the highly populated urns as the condensed phase. 
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Figure 3: The quantities e(t), and £ a (t) are plotted against t in a double logarithmic plot. The dashed lines are the 

analytically known equilibrium values e eq (Ti) and e eq (Tf) at the initial and final temperatures, respectively. In the inset the 
variances E(£), E £T (t) are shown together with the sum of the last two. 


t — 10 3 . e a {t) decays to zero proportionally to (In t)/y/t, as can be checked using the scaling form (05l) . 
showing that the internal energy is an equilibrating quantity and that, for large time, the contribution of e a 
can be neglected in e. 

Notice that the property lim^oo £^(t) = e eq (Tf) is an independent confirmation of the reliability of the 
method proposed above to separate fast and slow degrees. 

We have also computed the variance E(t) = (V. 2 ) — (' H ) 2 of the energy and of its gaseous and liquid 
parts, E^,(i) = (W. 2 )^ — (’%)(/, and E CT (f) = ('H 2 ) cr — (%)£, respectively. These quantities are shown in the 
inset of Fig. [3] Here one sees that for sufficiently long times (t > 500), when the separation of degrees 
becomes effective, E(f) = E^,(£) + E a (t). This is a nice indication that the i/)’s and the u’s are statistically 
independent variables. 


5.3. Autocorrelation function 

As an example of two-time quantity we consider the autocorrelation function. The number-number 
correlation 


C(t, t w ) = {n(t)ri(t w )) - p 2 


(28) 


which does not depend on i due to homogeneity, was analytically worked out in 22]. However this quantity 
is not well suited to the purposes of our paper because it is entirely dominated by the contribution of the 
highly populated urns forming the condensate. For this reason, instead of C , we will consider 


C(t,t w ) = {Sri(t)Sri(t w )}, 


(29) 


where Sr.i(t) = rj(£ + 1) — rj(f), is the variation of the number of marbles in the i -th urn. This quantity is 
related to C by 


C(£j t w ) — C(t + 1, t w + 1) — C{t + 1, t w ) — C(t , t w + 1) + C[t , t w ). 


(30) 
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Figure 4: The autocorrelation C(t,t w ) (continuous green curves) is plotted against t — t w on a semi-log scale, for three different 
waiting times (lower curves correspond to larger t w ). The long-dashed blue curve is C^(r) and the dot-dashed red one is 
t w ). In the upper inset the data collapse of the same curves of the left panel is obtained by plotting t^/ 2 C(t,t w ) against 
(t — t w )/t w . In the lower inset —C^,(t,t w ) is plotted for different choices of t w (see key). The cyan dots denote — C e q(T). 


The definition m is inspired by the autocorrelation encountered in the study of ferromagnetic systems. 
Typical values taken by the Sr's are 8ri(t) = 0,±1 (any value between — M and M is in principle possible, 
but with a very small probability) for all the urns, and this makes the difference with the definition (1281) 
where the largest contribution is provided by large values of r,; from the urns in the liquid phase. 

Let us now discuss the splitting m of this autocorrelation. The gaseous part has very different 
features with respect to the liquid one C a . If an urn is in the gas phase, the number of marbles contained 
in it fluctuates rapidly around the average gas occupation number p c = £(/3 — 1 )/£(/?) — 1 — 0.019 (where £ 
is the Riemann’s zeta-function), hence is expected to converge quite rapidly (over a microscopic typical 
time of the equilibrium state) to zero. This decay occurs from the negative side, because if by chance an 
urn with an average number of balls receives an extra marble, it will most probably release it in the future 
in order to keep its occupation number close to the mean. Furthermore, on general grounds, is expected 
to depend only on the time difference r. 

Conversely, in the effort to build the condensate, one urn in the liquid phase have the tendency to 
increase its number of marbles. This occurs up to a long time (of order t w ) when the competition with 
another prevailing urn will become effective and it will be progressively depleted. As a consequence, C a is 
expected to show a pronounced maximum at a time monotonously increasing with t w and then to decrease 
to zero. 

All these features are clearly displayed in the main part of Fig. []J where the autocorrelation is plotted 
against t — t w for different values of t w . Let us start by considering the solid green lines, which correspond 
to the whole correlation C(t,t w ). 

The curves relative to different values of t w collapse for short times, up to t — t w ~ 2. In this regime 
C(t,t w ) is negative. What we are observing in this region is clearly the evolution of the fast part C^{t). 

From t — t w ~ 2 onward the curves grow to a positive maximum, which is delayed and suppressed as t w 
is increased, and then decrease to zero. This is the behavior expected qualitatively for C a . Furthermore, in 
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this regime the curves for different t w can be collapsed by plotting tw 2 C(t,t w ), against t/t w , as it can be 
seen in the inset of Fig. 01 This means that the scaling form 


C(t,t w ) 



(31) 


is obeyed, which is of the expected type m, with the behavior m of £ previously found and 6 = 3, yet 
another consequence of the decrease GH) of urns in the liquid phase. 

Let us mention that the form m is in agreement with the known behavior 


C(t, t w ) = (p- p c ) tl /2 $(t/t w ), 


(32) 


(where $(cc) is a scaling function) of C(t,t w ) in the aging regime. Indeed, by plugging this expression into 
Eq. (1301) and expanding for large t w the various correlations on the r.h.s. around the value (t,t w ) of their 
arguments, one finds Eq. (I5T1) with c = — $'(t/t w ), where $'(a;) = d$>(x)/dx. 

In order to check the decomposition (fl3l) we have separated the urns into the gaseous and the liquid 
phase as discussed previously, and we have separately computed and C a by correlating among them 
degrees that pertain to one or the other set both at t w and t. These quantities are also plotted in Fig. 01 
where it is seen that they reproduce the whole quantity C in the short-time (C^) and in the long-time regime 
(C a ), respectively. For C a we have also checked that the scaling (1311) is obeyed. 

In order to better compare the behavior of the autocorrelation in the present urn model with that of a 
ferromagnet we have represented it pictorially in the right panel of Fig. [T] using the same line styles and 
colors in the two cases. The analysis of the two panels shows that a similar scenario occurs, apart of course 
from the different functional form of the various correlations. 

According to the general paradigm discussed in Sections [1] 01 01 C^(t) should correspond to the corre¬ 
lation function of the target equilibrium state. In the lower inset of Fig. 0 ] we plot this quantity computed 
in the aging regime by separating the degrees of freedom for different choices of t w . The data do not show 
any systematic dependence on t w , confirming that our definition of the set {ip} has the expected stationary 
property. In order to compare with the equilibrium autocorrelation we have computed the latter nu¬ 
merically: An efficient way to equilibrate the system at T <T C is to break the urn symmetry by starting 
from an initial configuration with all the marbles in one urn. This avoids the aging regime caused by the 
competition among urns and allows the system to quickly reach the equilibrium state. The quantity C eq (r) 
measured in this way is shown in the lower inset of Fig. 0 ] with a heavy-dashed cyan curve. As expected, 
C$(t) and C eq {r) coincide within numerical errors. The equilibrium autocorrelation could also be computed 


analytically starting from the known 2.2] expression of C 


C e q (t - ) — A e q,(/?y)r 


— (/3—3)/2 


(33) 


where A eq ((3) is a temperature-dependent constant. From this it is trivial to derive the expression 

C eq {r) * - (/3 ~ 1) 4 (/3 ~ 3) A e ,( / 3 / )r-^+ 1 )/ 2 (34) 

for C eq (t). We have checked that our numerical determination shown in the lower inset of Fig. 0 ] coincides 
within errors with Eq. (0H). 

Summarizing, our numerical data for the zeta urn model fit nicely in the general framework of slow 
relaxation discussed in Sections [T] El E] if the degrees of freedom of the system - the occupation numbers r* 
- are divided into two sets {ip} and {a} according to the recipe described in Section [foil 

Before closing this section, we briefly comment on the robustness of the overall picture with respect to 
different definitions - but same in spirit - of the fast and slow degrees. We tried to address this issue by 
changing the definition of ip's and a’s and, next to the technique of Section PTTTl we have separated slow and 
fast degrees using the running average procedure Eqs. ( 0161 ) of Section 0 ] More precisely, considering the 
autocorrelation function C we have defined the slow part as 


C a {t,t w ) = {Ar i(t) Ari(t ul )) 
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(35) 





where A rfft) = (A t)~ 1 [ri(t + At) — rj(f)] is the variation of number of marbles (normalized by At) in the 
i-th urn in a time At such that - as discussed in Section [3] - C^(t = At) ~ 0. Inspection of Fig. [I] suggests 
At ~ 10 4- 50 as reasonable values. By computing the quantity of Eq. (l35l) we found that, for sufficiently 
large values of t — t w the curve obtained superimposes on the determination of C a obtained previously with 
the technique of Section [5711 This results confirms the robustness of the degree-separation scenario. 


6. Conclusions 


In this paper we have discussed the paradigm of a simple class of aging systems where the dynamical 
properties can be accounted for by the coexistence of slow and fast degrees with different scaling properties. 

We have studied - in particular - the kinetics of the zeta urn model following a quench from high 
temperature to a final one below the condensation temperature. This model is a useful playground where 
fast and slow degrees have an intuitive interpretation in terms of the occupation numbers of different urns, 
allowing a distinction between a gaseous and a liquid phase, corresponding to the two sets {ip} and {cf} of 
fast and slow modes. 

We have considered the behavior of one-time quantities, as the energy, and of two-time ones, specifically 
the autocorrelation function. Due to the progressive reduction of the measure of the aging set {<r}, one-time 
quantities rapidly attain the characteristic value of the final equilibrium state. Non-equilibrium features are 
better observed in two-time quantities where the gaseous and the liquid component can be evidenced by 
increasing the two times along particular directions in the (t,t w ) plane, namely focusing on the short-time 
or on the aging regime, very similarly to what is known for ferromagnetic or glassy spin systems. 

The concept of an effective temperature is deeply related to the topic considered in the present Article. 
This quantity, which was originally introduced in the context of mean-field glass models [23-27], can be 


obtained from the relation between an autocorrelation function and the correspondent response function. 
For systems - like the one considered here - where two sets of degrees can be identified, one would expect 
two different temperatures associated to them. Accordingly, while the set {ip} is equilibrated at the bath 
temperature, a different effective temperature would characterize the slow part {<r}. An effective temperature 
for the zeta urn model based on the correlation (l28l) has already been investigated [13] : The different insight 
that could be provided by starting from our new definition (EH1) instead remains an open issue for future 
research. 
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